Effect of Dimple Potential on Ultraslow Light in a Bose-Einstein Condensate 
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We investigate the propagation of ultraslow optical pulse in atomic Bose-Einstein condensate 
in a harmonic trap decorated with a dimple potential. The role of dimple potential on the group 
velocity and time delay is studied. Since we consider the interatomic scattering interactions nonlinear 
Schrdinger equation or Gross-Pitaevskii equation is used in order to get the density profile of the 
atomic system. We find large group delays of order 1 msec in an atomic Bose-Einstein condensate 
in a harmonic trap with a deep dimple potential. 
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I. INTRODUCTION 

The impressive demonstration of ultraslow light prop- 
agation through an atomic Bose-Einstein condensate 
(BEC) l], utilizing electromagnetically induced trans- 
parency (EIT) 0], has promised an appealing applica- 
tion as a quantum memory Q. Beside the ultracold 
atoms, slow optical pulses have been observed in var- 
ious media like hot rubidium vapor 4]. In BECs, ul- 
traslow optical pulse is delayed by the order of few 
microseconds [l|. After, large controllable time delays for 
such broadband pulses were proposed Nevertheless 
large group delay of light was observed in hot rubidium 
gases [J . Quite recently, controlling the group velocity 
has been discussed @. It has been theoretically and ex- 
perimentally shown that this delay time can be increased 
by increasing the atomic density Q- Moreover, dimple 
potential can be used for increasing the time delay since 
the density of Bose-Einstein condensate can be increased 
by means of the dimple potentilas|8[. Phase space den- 
sity can be enhanced by an arbitrary factor by using a 
small dimple at the equilibrium point of the harmonic 
trapping potential Q. Recently, such potentials are also 
proposed for efficient loading and fast evaporative cooling 
to produce large BECs [Tfjj . 

If the atom-atom interactions can not be neglected, the 
structure of ground state of BEC is described by Gross- 
Pitaevski equation [ll[. If the scattering length a s is 
much less than the mean interparticle spacing, Gross- 
Pitaevski equation describes the zero-temperature prop- 
erties of the non-uniform Bose gas. 

In this paper, we investigate the one dimensional prop- 
agation of ultraslow optical pulse in an atomic Bosc- 
Einstcin condensate in a harmonic trap decorated with 
a dimple potential which is located at the center of har- 
monic potential. We study the role of the dimple poten- 
tial on the group velocity and the time delay. We model 
the dimple potential by a Gaussian function which has 
a narrow width value. Our calculations show that it is 
possible to observe large group delays up to 1 msec for 
an optical pulse in an atomic Bose-Einstein condensate 
which is trapped in harmonic potential decorated with a 
deep dimple. The paper is organized as follows: First of 



all the calculation of density profile of an atomic BEC 
system is briefly reviewed. After that, propagation of ul- 
traslow light under EIT system is discussed. We present 
our calculations and discuss the results in result section. 
Finally, we conclude in the last section. 



II. DENSITY PROFILE OF A TRAPPED BOSE 
GAS WITH A DIMPLE POTENTIAL 

The density profile of a Bose-Einstein condensate in 
an external potential can be approximated very well by 
means of Thomas-Fermi approximation. The density of 
ground state is given by the absolute square of the ground 
state wave funvtion: n(r) — \ip Q (r)\ 2 . Thomas- Fermi 
approximation (TEA) is often a good approximation to 
assume the total ultracold atomic density is in fact con- 
stant during the weak light field propagation. On the 
issue of the dimple, the validity of TEA depends on the 
length scale of the dimple. If it is much larger than the 
healing length, then TEA with the trap+dimple potential 
would work fine. The healing length or coherence length 
of the BEC is given by £ = (1/8-Knas) 1 / 2 [ijj]. Here n is 
the density of the ultracold atomic medium and we can 
take it as the peak density (n = p(0)) of the atomic sys- 
tem. We consider range of parameters in this work within 
the range of validity of TEA. One can approximate to- 
tal density profile of an ultracold atomic system for ID 
(ground state density) by [13] p(z) = [(/z — V(z))/Uq\. 
Here Uq = 4Trh 2 a s /m where m is atomic mass and a s is 
the atomic s-wave scattering length, /i is the chemical po- 
tential and can be evaluated by using Thomas-Fermi ap- 
proximation. The chemical potential is determined from 
N = J dzp(z). We consider the dimple potential which 
is modelled by a Gaussian function. We represent one 
dimensional harmonic potential with a dimple as 
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where ui z is the trap frequency of the harmonic trap in 
the z direction, and Vq > shows the strength (depth) 



of the dimple trap (V^ 
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It is possible to increase Vq from to 1500fi.w z [l3|. We 
use large values for the width of Gaussian function l z (in 
G(z) — Vo exp (— z 2 /2l 2 )) compared to the extension of 
BEC (see Sec IV.) in order to get an appropriate mod- 
eling for the narrow dimple. If we apply a deep dimple 
to the atomic condensate in a harmonic trap we will get 
the total number of atoms analytically by: 



N 



2(m+Vq) 
(^i+Vo/lf) 1 

az — 

U 

x {fi- ^muj 2 z z 2 + V e~ i: ^ )2 }. (2) 



By taking the integral in Eq.@, we find an analytical 
expression for total number of atoms in terms of chemical 
potential, interaction term, trap frequency and strength 
of the dimple potential: 



N = 



: {2V0« + Vrj)(3/iVb 



3V2Uo(V + muj 2 l 2 ) 3 / 2 
+ mu 2 l 2 z (2fi-Vo))+3Vo(V + mw 2 z l 2 z ) 

x y/(V + mufflirErfL L, M \ \ n )}. (3) 



Vq + muj 2 l 2 ' 



Here erf(x) is the error function. The error function 
can be expanded in terms of x where x = yj Vo+tnu kl?' 
erf(x) = -L.[2x - ^ + ^ - + 0[x} 9 }. We insert 

this expanded term into Eq.Q and solve this equation 
numerically in order to get chemical potential. Doing this 
we find the value of the chemical potential p = 1.0800 x 
10~ 12 eV for m = 23 amu ( 23 Na), ui = 200 Hz and V = 
100huj z = 6.5821 10"13 eV for JV = lx 10 6 . 



III. ULTRASLOW LIGHT UNDER EIT SCHEME 

We consider an EIT model for a gas of TV three-level 
atoms interacting with two laser beams in A configura- 
tion. The upper level is coupled to the lower levels via 
a strong drive field with frequency w c and a weak probe 
field of frequency lu p . At resonance the absorption of the 
probe field can be neglected. Weak probe beam propa- 
gates along the condensate axis in the z direction. Prop- 
agation of the ultraslow wave packet in one dimensional 
inhomogeneous atomic condensate can be described by 

MM 



dE 

dz 



■a{z)E- 



v g (z) dt 



dE , . ,d 2 E 
— +*6 2 (z) — 



0, (4) 



where a{z) is the pulse attenuation factor, v g (z) is the 
group velocity, and 6 2 (z) is the group velocity disper- 
sion. The third order dispersion is found to be much 
smaller and neglected [Hj]. EIT susceptibility [l6| of Bose- 
Einstein condensate of atomic density p is expressed as 



X = PXi with 

H 2 



Xi 



i(-iA + r 2 /2) 



e h [(T 2 /2 - iA)(r 3 /2 - iA) + tt 2 /4] ' 



(5) 



where A = to — loq is the detuning of the probe field 
frequency oj from the atomic resonance ujq. In Eq. ([5]), 
f2 c is the Rabi frequency of the control field and /i is 
the dipole matrix element for the probe transition. T2 
and T3 denote the dephasing rates of the atomic levels. 
The significant position dependent group velocity for the 
optical pulse propagation can be calculated from the sus- 
ceptibility using the relation[14j 



(6) 



Here we take A = 589nm, and ujq = 2irc/\. In Eq.© 
group velocity depends on atomic density. As mentioned 
in the introduction, the atomic density can be increased 
by using a dimple potential Q. An optical ultraslow 
pulse propagates through the ultracold medium without 
absorption due to the small imaginary part of the EIT 
susceptibility at resonance. 



IV. RESULTS AND DISCUSSIONS 

We consider a gas of N = 1 x 10 6 23 Na atoms with 
T 3 = O.57, 7 = 6 MHz, T 2 = 6 x 10 3 Hz, and n c = 
O.57. We take accessible experimentally parameters such 
as u r = 350 Hz and uj z — 100 Hz, so that peak density 
p a = 1.56 x 10 20 1/m 3 . The Eq.® is solved numerically 
in order to find the chemical potential. We assume that 
all atoms are loaded into the harmonic potential with a 
dimple therefore effective length of the atomic medium 
becomes at the order 6/im. The density of the condensate 
is mainly controlled by the dimple potential for extremely 
deep dimple. We present the change of p, as a function 
of Vq in Fig.([T]). As the strength of the dimple potential 
increases, chemical potential becomes larger. 

When light enters the condensate, its group speed ex- 
hibits a dramatic slow down. Here we consider resonant 
probe pulse with A = 0. Within the condensate region, 
at zero temperatures, the group velocity remains approx- 
imately at the same ultraslow value. Light rapidly accel- 
erates to high speeds when it leaves the condensate at 
the interface to thermal part. At extremely low temper- 
atures, po saturates to Thomas-Fermi density. At zero 
temperature, vgo can be calculated by po which was con- 
verted from Thomas-Fermi approximation. We can ig- 
nore the spatial variations therefore density of the atomic 
system can be taken as p(z = 0) = po- in other words, 
at low temperatures group velocity is determined by po 
and group velocity decreases rapidly with the increas- 
ing strength of the dimple for small dimple strengths as 
seen in Fig. ((2]). However, group velocity decreases slowly 
at deep dimple. One can interpret the decrease in the 
rate of slow down as follows: As Vq becomes larger, the 
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FIG. 1: The chemical potential fi vs strength of dimple po- 
tential for 23 Na the Bose-Einstein condensate of TV = 1 x 10 6 
atoms. The dots show the chemical potential which corre- 
spond to the strength of the dimple potential. Solid line rep- 
resents the linear fitting. Dimple potential (Vb) is scaled by 
Tiu) z . The parameters used are M = 23 amu, a s — 2.75 nm. 
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dimple part predominates over the harmonic part of the 
potential. 

The group velocity is almost constant within the ultra 
cold atomic medium. In this case we ignore the small 
contributions of modal and waveguide dispersions and 
determine the group velocity, same for both fractions, by 
assuming a constant peak density of the condensate in 
the material dispersion relation. 

So the delay time can be calculated approximately by 
the time delay formula tjj = L z /v g . Here tjj is the time 
delay of the ultraslow pulse, v g is the group velocity in 
the z direction and L z is the axial length of the conden- 
sate. Here axial length of the condensate (L z ) can be 
taken as: L z — 2R where R is the Thomas-Fermi axial 
radius which is given by y / 2(/i + Vo)/muj^. Therefore we 
find large group delays at the order 1 msec in an atomic 
Bose-Einstein condensate in a harmonic trap with a deep 
dimple potential in which the strength of the deep dimple 
potential is Vn = 1500?ia; z . 

V. CONCLUSION 

We have explored the propagation of ultra slow light 
through a Bose-Einstein condensate in a harmonic trap 
with a dimple potential. We have investigated the ef- 
fect of the dimple potential on the group velocity and 
the time delay. As the strength of the dimple potential 
increases group velocity becomes smaller. However, at a 
critical value of the dimple group velocity becomes ap- 
proximately constant. As a result, time delay can be in- 
creased by means of an atomic Bose-Einstein condensate 
in a harmonic trap with a dimple potential. 



FIG. 2: Group velocity vs strength of dimple potential for 
A = 0, propagating through a 23 Na Bose-Einstein condensate 



under EIT scheme for TV 



1 x 10 atoms. The solid line 



shows the peak value of the group velocity vs strength of the 
dimple potential. Dimple potential (Vb) is scaled by huj z . 
The parameters are the same those of Fig[T] 
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